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Effects  of  Shear  and  Strain 

1.  Won-Uniformities. 

An  understanding  of  the  response  of  a  pre-mixed  flame  to  non-uniformities 
in  the  flow  is  important  in  many  technological  situations.  To  sustain  a  flame 
in  a  high-velocity  stream  the  turbine  engineer  must  provide  anchors,  which 
generate  strong  shear.  The  designer  of  an  internal  combustion  engine  is 
concerned  with  burning  rate  in  the  swirling  flow  above  the  piston.  Turbulence 
is  all  pervasive  and  then  the  flame  is  subject  to  highly  unsteady  shear  and 
strain.  These  situations  are  extremely  complicated,  of  course,  and  it  is 
unlikely  that  mathematical  analysis  will  ever  provide  detailed  descriptions; 
those  must  be  left  to  empirical  studies  augmented  by  massive  numerical  computa¬ 
tions.  Nevertheless,  analysis  of  the  response  of  a  flame  to  a  simple  shear, 
for  example,  can  provide  useful  insight  into  the  interaction  mechanism  in 
more  complex  situations.  Indeed  the  answers  to  such  elementary  questions  can 
be  used  as  a  much-needed  guide  in  studies  of  more  complicated  ones. 

Moreover  there  are  simple  circumstances  in  which  such  analysis  has  direct 
significance.  A  burner  flame  experiences  shear  in  the  neighborhood  of  the  rim, 
ahd  its  quenching  depends  on  the  local  character  of  that  shear.  A  flame 
immersed  in  a  laminar  boundary  layer  experiences  both  shear,  due  to  velocity 
variations  across  the  layer,  and  strain,  due  to  streamwise  variations,  and  its 
quenching  will  depend  on  their  local  values. 

For  a  constant-density  flow  the  relative  motion  near  a  point  is  the 
superposition  of  a  simple  shear  and  two  simple  strains  (cf.  Batchelor  19&7, 
p.  79).  These  two  elementary  flows,  simple  shear  and  simple  strain,  are 
therefore  basic  to  any  discussion  of  the  effect  of  non-uniformities  on  flames; 
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the  present  chapter  is  for  the  most  part  devoted  to  them.  Only  premixed 
flames  are  discussed  and,  except  at  the  end  (Sec.  5),  the  formulation  adopted 
is  that  of  near-equidiffusional  flames  as  simplified  in  Sec.  VIII. 6. 

The  fundamental  question  is  whether  or  not  a  flame  can  be  extinguished 
by  a  non-uniform  flow  per  se  .  In  Sec.  IX. 5  it  was  shown  that  the  proximity 
of  a  cold  wall  can  quench  an  equidiffusional  flame,  the  mechanism  being  heat 
loss  from  the  combustion  field  rather  than  the  geometrical  constraint  imposed 
on  it  by  the  wall.  If  flow  non-uniformities  can  induce  a  similar  heat  flux, 
quenching  will  occur  with  the  cold  unburnt  gas  replacing  the  wall  as  heat 
sink.  But  even  if  they  cannot,  we  may  still  expect  unequal  diffusion  of  heat 
and  reactant  (i.e.  X  5*  o)  to  provide  other  circumstances  for  which  there  is 
quenching,  since  that  is  the  case  for  the  flame  tips  of  Sec.  IX.j*.  When  X 
was  sufficiently  negative,  the  presence  of  a  boundary  quenched  plane  tips  and 
flame  curvature  quenched  axisymmetric,  w'thout  the  help  of  heat  loss. 

2.  Response  to  Simple  Shear. 

Consider  a  plane  parallel  flow  which  is  uniform  for  y  <  0  but 
linearly  sheared  for  y  <0,  i.e. 

(1)  V  =  (Uf ,0) 

where 

]1  for  y  >  0, 

(2)  f(y)  =\ 

(^l-ioy  for  y  <  0. 

In  the  limit  U  »  the  equations  with  which  we  have  to  deal  are 

(3)  f9(T,h)/3x=  kQ(T,h  +  XT)  for  -«  <  y  <  F(x), 
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(4)  T  =  T#  =  H1,  f  3h/3x  =  LQ(h)  for  F( x)  <  y  <  », 

these  being  the  modifications  of  equations  (lX.31,32)  which  take  account  of 
the  flow  (2).  With  them  go  the  Jump  conditions  (IX. 34),  together  with  the 
boundary  conditions  (IX. 36)  and 

(5)  T  -*•  T- ,  h  0  as  y  • 

The  uniform  part  of  the  flow  can  support  a  steady  plane  deflagration 
wave  inclined  at  an  angle  such  that  the  normal  component  of  the  gas  speed 
equals  the  adiabatic  flame  speed,  and  this  is  assumed  to  be  the  form  of  the 
combustion  field  for  y  large  and  positive  (Fig.  l).  In  other  words,  the 
initial  conditions  (IX. 29, 30)  still  hold,  and  indeed  the  first  for 
_®  <  y  <  f(0).  As  the  flame  approaches  the  x~axis  it  is  influenced  by  the 
shear  as  soon  as  its  preheat  zone  penetrates  the  lower  half-plane.  The 
increased  gas  speed  there  may  be  expected  to  deflect  the  flame,  so  that  its 
slope  decreases. 

The  above  formulation  is  due  to  Buckmaster  ( 1979a),  who  reaches  the 
following  conclusions.  For  some  values  of  the  shear  gradient  u  and 
Lewis-number  parameter  X,  the  slope  actually  decreases  to  zero.  At  that 
point  the  flame  speed  is  zero  and,  according  to  the  hypothesis  of  Sec.  IX. 4, 
quenching  occurs.  Invariably  the  quenching  point  lies  in  the  uniform  region 
although  the  preheat  zone  has,  of  course,  penetrated  the  shear  region.  For 
the  remaining  values  of  <o  and  X  the  flame  is  never  quenched  even  though  it 
may  fail  to  penetrate  the  shear  region.  A  rationalization  for  not  being 
quenched  in  the  shear  region  is  that  the  non-uniformity,  as  measured  by  f'/f, 
is  a  maximum  at  y  =  0  so  that,  if  the  flame  manages  to  cross  the  boundary, 
it  can  .survive  the  ever-weaker  non-uniformity  beyond.  Certainly  there  is  a 
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solution  of  the  governing  equations  as  y  ■+■  -«  corresponding  to  a  locally 
plane  flame  propagating  with  unit  speed  (cf.  Sec.  9.5). 

Clear  evidence  of  the  dichotomy  is  shown  for  m  =  5  by  the  numerical 
results  plotted  in  Fig.  2.  For  X  =  0  or  5  the  flame  penetrates  the  shear  and 
shows  no  sign  of  quenching;  indeed  no  solution  showed  the  flame  speed 
approaching  zero  in  y  <  0.  But  for  X  =  10  the  flame  speed  falls  to  zero 
at  the  point  Q  outside  the  shear  region.  The  figure  also  shows  that  the 
shear  has  more  effect  on  the  flame  as  X  increases,  in  contrast  to  flame 
tips  (Fig.  IX. 6)  where  the  boundary,  i.e.  symmetry  line,  has  less.  Moreover, 
the  present  trend  does  not  agree  with  that  for  slowly  varying  flames  in  Sec. 

The  effect  of  increasing  to  can  be  seen  from  the  limit  w  -*  ».  There  i 
then  a  fully  developed  combustion  field,  i.e.  an  asymptotic  solution  as 
X  -*■  »  ,  in  which  the  flame  lies  at  a  definite  level  in  y  >  0  and  T,h 
are  independent  of  x  •  In  fact,  in  the  shear_ region  away  from  the  boundary 
convection  dominates  for  all  x  »  so  that 

(6)  3T/8x  =  3h/3x  =  0  for  y  <  0, 

i.e.  there  is  no  dependence  on  x  anywhere.  The  problem  is  therefore  reduce 
to  one  in  y  >  0  only,  under  the  boundary  conditions 

(7)  T  *  Tx,  h  =  0  on  y  =  0. 

In  the  limit  x  *►  05  the  initial  conditions  may  be  ignored  and  the  system  then 
has  the  solution 

(8)  T  =  T1  +  Y1y/F(oo),  h  =  -  XY^/FCoo)  for  0  <  y  <  F(»), 

#  =  H^,  h  =  -XY1  for  F(co)  <  y  <  «, 


(9) 


T  =  T, 
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vhere 

(10)  P(»)  =  exp(XY1/2T#) 

is  the  asymptotic  height  of  the  flame.  The  expression  (9b)  violates  the 
boundary  condition  (IX.  36)  for  h,  showing  that  the  limits  y  -*•  »  and  X  00 
do  not  commute.  There  is  an  adjustment  at  large  values  of  y,  which  depends 
on  the  history  of  the  flame;  but  that  is  of  no  importance  for  our  purposes. 

F(")  is  large  even  for  moderate  values  of  X.  On  the  other  hand, 
inspection  of  Fig.  2  shows  that  flames  are  first  affected  by  the  shear  when 
F  is  about  2.  The  implication  when  m  =  »  is  that,  for  sufficiently  large 
X  >  0,  there  is  a  minimum  in  F,  i.e.  the  flame  speed  must  have  fallen  to  zero 
at  some  finite  value  of  y.  Numerical  solutions  of  the  limit  problem  (Fig.  3) 
point  to  X  =  0  as  the  critical  value:  for  X  <  0  the  function  F  decreases 
monotonically  to  its  asymptotic  value  (10),  whereas  for  X  >  0  it  has  a 
minimum  (which  we  have  taken  to  imply  quenching).  We  conclude  that  a 
sufficiently  strong  shear  flow  will  prevent  a  flame  from  reaching  the  boundary 
y  =  0  and  actually  quenches  it  for  X  >  1. 

Eventually  the  limit  equations  break  down  and  a  new  coordinate  y/w  is 
needed,  corresponding  to  significant  diffusion  in  the  shear  region.  Bucknaster 
(1979a)  argues  that  F  increases  again  when  X  <  0  is  sufficiently  large, 
which  implies  quenching  for  X  sufficiently  smaller  than  one  also.  However, 
no  numerical  results  exhibiting  such  quenching  have  yet  been  obtained. 

3.  Response  to  Simple  Strain. 

In  contrast  to  simple  shear,  the  variation  of  the  essential  velocity 
component  in  simple  strain  is  in  its  own  direction  rather  than  at  right  angles. 
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The  distinguished  limit  of  Sec.  VIII. 6,  in  which  the  non-dimensional  gas 
velocity  and  the  length  scale  on  which  it  varies  become  unbounded  but 
their  ratio  remains  fixed,  leads  to  the  governing  equations  (VIII.70).  The 
unit  of  the  velocity  gradient  represented  by  €  is  that  of  the  flame 
speed  divided  by  the  thickness  under  adiabatic  conditions.  As  we  shall  see, 
a  velocity  gradient  that  is  small  on  this  scale  has  little  effect  on  the 
flame,  whereas  a  large  enough  gradient  extinguishes  it. 

Equations  (VIII.70)  hold  ahead  of  the  flame,  i.e.  for  y  >  F(x).  Behind 
the  flame  we  have 


(ll)  T  =  T#  =  H  ,  h  =  h#  for  0  <  y  <  F(x) 

if  (a)  the  wall  is  insulating  or  (b)  the  flame  is  located  in  equal  opposing 
jets  ( Fig .  U)  so  that 


(12)  3h/3y  =0  on  y  =  0. 


(The  condition  ensures  that  wall  quenching  does  not  intrude.)  For  case  (a) 
the  flame  will  in  reality  be  imbedded  in  a  boundary  layer  (unless  the  Fran:;;], 
number  is  very  small),  so  that  the  model  is  more  directly  relevant  to  case  (b). 
We  must  add  the  boundary  conditions 


(13) 


T  +  h  +  0  as  y  -> 


and  th  •  jump  conditions  (IX. 3^). 

Since  x  only  appears  in  F  there  is  a  solution  with  the  flame  lying 
parallel  to  y  -  0,  the  constant  value  of  7  being  given  implicitly  by 


0«0 


-dc 


e 


'V  °d  erfc  d  expt-^-T^-2 -  -  -  d^) ] , 


2ll^  7r  ^erfc  d 


J 
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where 


B  =  (6/2)  ,  d  =  8F. 


The  solution  ahead  of  the  flame  is 


T  =  T^  +  erfc(8y)/erfc  d 


(17)  h  =  h* 


~erf c^~d~~  +  ~h — ^~2~  ^ye  ^  Y  erfc  d  "  de_d  erfc(By)l 
it  2erf  c  d 


for  y  >  F. 


where 


2 

(18)  h*  =  2H^n(2g  e"d  /  Arfc  d). 

This  determination  of  h*  also  completes  the  description  (ll)  "behind  the 
flame . 

Of  greatest  interest  i$  the  variation  in  stand-off  distance  F0  with 

the  straining  rate  £  (Fig.  5).  There  are  two  kinds  of  response,  depending 

2 

on  the  value  of  X.  For  X  <  4^,  it  is  single-valued,  with  Ffl  decreasing 
monotonically  to  zero  at  £  =  € where 

(19)  £  =  it  exp(-X/2H!r) . 

c  J* 

For  £  >£,,  there  is.no  solution.  For  X  >  tK^,  the  response  is  double-valued 
for  some  values  of  £  ;  and  there  is  no  solution  beyond  a  value  6  ,  at  which 

x=Fe, 

F ViS'non-zero.  Corresponding  values  of  £g  and  F£  for  various  values  of 
X  are  listed  in  Table  1  for  the  same  parameter  value  H  =  1.2  as  in  Fig.  5. 


A 


Table  1 


X 

5.76 

6.03 

6.48 

7.03 

7.52 

8.05 

9.05 

10.02 

11.02 

12.02 

l4. 08 

16.00 

18.  CO 

F 

e 

0 

.23 

.62 

1.11 

1.57 

2.07 

3.05 

4.03 

5.05 

6.10 

8.29 

10.36 

12.55 

fe 

.212 

.194 

.170 

.146 

.130 

.116 

.096 

.082 

.072 

.064 

.052 

.044 

0*5  P 

The  physical  picture  is  now  clear.  Suppose  the  strain  is  increased 
from  a  small  value  sufficiently  slowly  for  the  combustion  field  to  be  quasi¬ 
steady.  The  increase  causes  the  flame  to  move  closer  to  the  stagnation  point. 
If  X  is  large  enough,  it  reaches  the  stagnation  point  for  the  critical  value 
6^ ,  and  is  extinguished.  Otherwise  extinction  occurs  for  €  ,  when  it  is 
still  a  distance  Fg  away.  In  the  latter  case  the  lover  branch  of  the 
response  (where  F^  increases  with  f  ) is  unstable,  as  has  been  shown  by 
Buckmaster  ( 1979b).  (His  treatment,  docs  not  rule  out  other  parts  of  the 
response  being  unstable  too,  so  that  the  extinction  picture  above  may  have 
to  be  modified. ) 

The  flame  speed  W  (Fig.  6)  i.e.  the  ga  speed  normal  to  the  flame  sheet, 

has  the  value  s'  p.  As  €  ->•  0,  so  that  F  ->  =»  ,  it  tends  to  the  adiabatic 

o 

value  1;  when  t  =  6  it  vanishes,  since  F  is  zero.  Like  it  is 

C  O  O 

2  2 
single-valued  for  X  <4ll^  and  partly  double-valued  for  X  >  4r  .  There  is  an 

2 

initial  increase  in  V/  for  X  <-2H^;  otherwise  it  decreases  monotonically 

until  C  =  f  or  (:  .  It  is  therefore  possible  for  the  flame  to  have  a 
c  e 

speed  greater  than  its  adiabatic  value,  even  though  it  is  being  stretched 
at  a  positive  rate. 

In  short,  there  is  a  limit  (depending  on  the  Lewis  number)  to  the  amount 
of  strain  that  a  premixed  flame  can  tolerate.  The  con  civ.,  i:u  is  consistent  vi 


the  experimental  fact  that  a  flame  held  at  the  front  of  a  bluff  body  can  be 
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blown  off  by  increasing  the  gas  speed  sufficiently.  Attachment  is  then 
usually  transferred  to  the  rear. 

k.  Response  to  More  General  Non-Uniformity. 

Reduction  to  a  parabolic  problem  occurs  for  general  non-uniform  flows 
in  which  the  velocity  is  everywhere  large  but  changes  proportionately  over 
correspondingly  large  distances  (Buckmaster  1979b).  An  example  is  the 
potential  motion  of  a  fast  stream  past  a  large  cylinder,  whose  front  stagnatio 
point  was  the  subject  of  the  last  section.  Such  flows,  when  plane,  are 
characterized  by  velocity  fields 

(20)  V  =  Ua(x/U,y/U); 

'letting  U  -*■  m  gives  a  distinguished  limit  in  which  ^  y  remains  finite. 

We  shall  see  that  the  effective  part  of  the  limiting  motion  is,  on  the  x,y-sca 
made  up  of  a  fast  uniform  flow  and  a  simple  strain,  both  of  which  vary  slowly. 
Shear  is  ineffective,  so  that  the  term  general  straining  motion  is  appropriate 
Take  curvilinear  coordinates  with  x  measuring  distance  along  one 
particular  streamline  and  y  distance  from  it.  Away  from  a  stagnation  point 
these  will  be  slowly  varying  cartesian  coordinates,  since  the  curvature  of 
streamlines  is  vanishingly  small.  The  angle  between  the  flame  sheet  and  an 
intersecting  streamline  is  also  small,  so  that 

(21)  X  =  x/U  and  y 

are  appropriate  variables  for  T  and  h.  Likewise  the  velocity  field  may 
be  approximated  by 
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(22) 


Y  =  (UqQ,  -  q£y). 


where  q^(x)  is  "the  speed  on  y  =  0  and  continuity  has  been  taken  into  account. 
Only  leading  terms  have  been  retained,  so  that  a  simple  shear  due  to  variation 
of  the  x-component  in  the  y-direction  is  ignored.  Uniformity  breaks  down  at 
large  distances  from  the  flame  sheet,  but  it  does  not  matter  there.  On  the 
x,y-scale  the  velocity  relative  to  that  at  any  point  is  a  simple  strain  with 
€  =  (Sec.  VIII. 6).  Together  with  the  fast  uniform  flow  (Uq^,0)  at  the 
point,  both  of  which  vary  with  x  >  it  makes  up  the  effective  part  of  the 
limiting  motion. 

In  the  limit,  only  the  approximation  (22)  enters  into  the  governing 
equations  (XIII.62,  63),  which  become. 


(23) 


(qo  £  -  *  Lo(T-  h  +  w)' 


^0  3y 

The  response  to  an  arbitrary  streamwise  velocity  q^(  )d  may  be  determined 
from  the  computer  as  easily  as  that  to  simple  shear  (Sec.  2).  Note  that 
equations  (VIII.67,  68),  which  were  used  in  Sec.  IX. U,  are  the  special  case 
=  1  when  the  strain  vanishes  everywhere.  The  simple  shear  flows  discussed 
in  Sec.  IX. 5  and  Sec.  2  are  not  of  the  type  considered  here  since  changes  of 
the  velociy  field  in  the  y-direction  are  large.  As  already  noted,  the  analysis 
breaks  down  at  a  stagnation  point,  so  that  equations  ( VIII. 70)  cannot  be  recoverec 
Nevertheless,  "he  stagnation-point  solution  in  Sec.  3  has  the  required  form  for 
large  x,  with  =  6  x>  and  hence  provides  initial  conditions  for  any 
q^  ~  £x  as  x  0  when  the  boundary  conditions  on  T  and  h  are  compatible. 
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Two  examples  will  be  considered,  namely 


(24) 

(25) 


qo(x)  = 


{€  sin  x  for  0  <  x  < 

'  €  (x~  2X3/3  +  x5/5)  for  X  > 


0. 


The  first  corresponds  to  potential  flow  around  the  cylinder  and  both  asymptote 


<r  x- 

If  therefore  the  same 

boundary  conditions  as  in  Sec.  3  are  imposed. 

namely 

(26) 

T  T  ,  h  -»■  0 

as  y 

(27) 

dh/Sy  =  0  on 

ii 

o 

S3 

then  the  results  obtained  there  provide  appropriate  initial  conditions  here. 

Of  course,  behind  the  flame  the  appropriate  solution  of  equations  (23)  has 
T  =  T#  =  Hr 

The  flows  (24)  and  (25)  are  quite  different  in  nature.  For  the  first 
the  stretch,  represented  by  the  Karlovitz  number  q^,  decreases  monotonically 
from  6  to  -€  (Fig.  7a);  but  for  the  second  it  decreases  from  €  to  a 
minimum  of  zero  at  x  =  ^  and  then  increases  indefinitely  (Fig.  7b).  These 
characteristics  are  particularly  relevant  to  an  examination  of  the  flame 
speed 

(28)  W=(qQF)', 

which  can  deduced  from  the  computations  by  numerical  differentiation. 

Fig.  7  shows  that  W  increases  monotonically  for  the  flow  (24);  but 
that  for  the  flow  (25)  it  increases  initially,  reaches  a  maximum  around  x  =  1» 
and  then  decreases.  It  is  hard  to  resist  the  conclusion  that  there  is  an 


inverse  correlation  between  the  changes  in  flame  speed  and  stretch,  at 
least  for  moderate  values  of  X.  The  correlation  is  not  precise  (as  we 
have  already  seen  in  Sec.  3):  in  Fig.  7b  the  flame  speed  achieves  its 
maximum  exactly  at  x  =  1  for  only  one  value  of  A  (which  appears  to  be 
zero).  But  the  general  trend  is  undeniable. 

The  computations  for  the  flow  (25)  show  that  the  decrease  in  W  from 
its  maximum  continues  until  zero  is  reached  where,  according  to  Sec.  IX. U, 
quenching  occurs.  The  results  therefore  suggest  the  general  conjecture  that 
a  flame  subjected  to  sufficiently  strong  and  increasing  strain  will  be 
extinguished.  Hitherto  a  sound  mathematical  basis  for  such  a  conjecture  has 
been  lacking. 

5.  The  Effects  of  Slowly  Varying  Shear  and  Strain. 

The  chapter  ends  with  an  examination  of  slowly  varying  flames  under  the 
influence  of  simple  shear  and  simple  strain.  For  strain  the  results  merge, 
as  o£  -*■  1  ±  0,  with  those  of  near-equidiffusional  flames,  as  A  -*+»;  but  for 
shear  they  do  not.  We  have  no  explanation  to  offer. 

For  simple  shear,  the  5 ,n-coordinates  of  Sec.  VIII. 3  are  now  taken  as 
shown  in  Fig.  8  .  Then  the  velocity  field  to  be  inserted  in  equation  (VIII.32) 
is 


(29) 

v  =  (Uf,0) 

where 

C  1  for  ri>  0 

(30) 

f(n)  = 

v  l-wri  for  n  <0. 

-13- 


If 

(31)  n  -  F(5) 

is  taken  to  be  the  locus  of  the  flame,  then  for  n  <0  we  have 

(32)  M  ■  p,U(l  -  uF)  |F'  I  (l  +  F,2)"\  v,  =  U(l-u)F)sgnF*(l  +  F*2)-^ 

n  ±  11  -i— l 

if  the  latter  is  measured  in  the  s-direction  shown.  Hence 

(33)  y±1  -VjMn  -  M^-v^  =  M2  sgn  F’(l  +  f'W/PjF’2 

and  the  basic  equation  ( VIII. 32)  becomes 

(3*0  F"  =  k(l-uF)F  3£ruM 

n 

where  k  is  the  constant  (IX.21). 

The  corresponding  formulas  for  n  >  0  are  obtained  on  setting  u  =  0, 
and  then  equation  (IX.20),  governing  the  shape  of  a  flame  tip  in  a  uniform 
flow,  is  recovered  (with  tan^  =  F')«  The  appropriate  solution  is 

(35)  F'  =  -tan  a 

corresponding  to  a  plane  flame,  althoxigh  a  curved  flame  (corresponding  to 
the  tip)  is  also  available  when  is  greater  than  1.  We  therefore  take 

(36)  F(0)  =  0,  F' (0)  =  -tan  a 

as  initial  conditions  for  the  governing  equation  (3M  in  n  <  0. 

Fig.  *?  shows  the  effect  of  the  shear.  For  J_<1  the  flame  is  bent 
towards  the  flow  direction  at  a  rate  that  increases  with  w  ,  but  is  never 
actually  quenched.  On  the  ether  hand,  the  flame  bends  away  from  the  flow 


-In¬ 


direction  for^_>  1,  again  at  a  rate  that  increases  with  to  .  These  results 

fa* 

are  consistent  with  those  obtained^flame  tips  in  a  uniform  flow  (Sec.  IX. 3), 
where  a  rounded  tip  was  possible  for  >1  but  not  for  <^_<  1.  However,  they 
do  not  merge  with  those  for  a  near-equidiffusional  flame  (Sec.  2),  where  the 
turning  ability  of  the  shear  increased  with  X.  (The  inconsistency  of  Secs.  2 
and  IX. 4  has  already  been  noted  in  the  later  section.) 

We  turn  now  to  simple  strain,  i.e.  stagnation-point  flow,  where  the  results 
do  merge  with  those  for  near-equidiffusional  flames.  The  comparison  involves 
the  so-called  flat  flames  here;  for  <  1  it  turns  out  that  such  flames  can 
tolerate  any  amount  of  strain,  but  for  ^  >1  there  is  a  limit. 

With  axes  as  shown  in  Fig.  10,  the  velocity  field  with  which  we  have 


to  deal  is 


Y  “  ^ (S,-n)- 


For  the  locus  ( 31 )  we  find 


(38)  Mn  =  p^UF'  +  F)/(l  +  f'2)^,  v_Ll  =  6U  -  FF')/(1  +  f'2)^ 


so  that 


YJ-1  ’  ~-L Mn  “  Mn  ~J_'Yll  =  p  (-L+p'S)^  +£F' } 

and  the  basic  equation  ( VIII. 32)  becomes 

(39)  (52  +  F2)F"  +  (1  +  f'2)(cF’  -  F)  =  k(F  +  £F')3£n  M  , 


where 


k  =  -2T#p![f/b. 


Solution  of  this  equation  requires  two  boundary  conditions  but  only 


one  is  obvious ,  namely  the  symmetry  condition 
(Ul)  F'(0)  =  0  . 

The  difficulty  reflects  the  elliptic  nature  of  the  problem  and,  for  determinacy, 
some  reference  must  be  made  to  the  far  field.  If,  for  example,  we  are  dealing 
with  a  fast  flow  past  a  correspondingly  large  cylinder  (i.e.  with  £  fixed) 
then  it  is  appropriate  to  look  for  a  constant  solution,  corresponding  to  a 
flame  that  varies  on  the  scale  of  the  cylinder.  Such  a  solution  is  given  by 

( k2 )  M^n  M2  =  -2p2£2/k  with  M  »  p,  £  F. 

n  n  ±  n  l 

Sivashinsky  (1977)*  who  considers  the  axisymmetric  counterpart  of  the 
present  problem,  propones  that  such  a  choice  is  correct  in  all  circumstances 
but  there  is  no  evidence  for  his  view.  Indeed,  completely  acceptable  solutions 
are  obtained  by  giving  F(0)  a  value  that  does  not  satisfy  equation  (b2) , 
see  Fig.  11.  A  constant  solution  is  just  a  separatrix  or  an  asymptote  in 
the  figure. 

Fig.  12  shows  that  there  is  always  a  unique  solution  of  equation  (1+2) 
for  b  >  0  but  that  the  condition 

(1+3)  £  <_  -T2/eb 

must  be  met  for  b  <  0,  and  then  there  are  two  solutions.  We  conclude  that 
there  is  always  a  flat  flame  for^<  1  but  that  there  is  a  limit  to  the 
straining  rate  for  "/,  >1.  To  reconcile  these  results  with  those  for  near 
equidiffusional  flames  note  that  here  €  is  measured  on  the  scale  G-'*'  so  that 


-l  6- 


the  limit  (^3)  tending  to  infinity  as  1  +  0  is  (b  -0)  is  compatible 
with  €  remaining  finite  there  as  X  -*■  00  (suggested  by  Table  l).  On  the 
other  hand,  -*•  ®  as  X  agrees  with  the  absence  of  a  limit  here  for 

any  X  >  1  (in  particular  «£.-  1  +  0). 

Fig.  12  could  have  been  introduced  after  the  general  relation  (VIII.l*3) 
to  suggest  that  extinction  will  occur  whenever 


m 


1  dA  <  21 

4  di  >  be 


for  b  <  0. 


However,  without  a  guarantee  that  such  voluminal  stretch  is  possible,  the 

speculation  would  have  been  idle.  Stagnation-point  flow  provides  the  first 

example,  and  so  far  the  only  one.  If  we  restrict  attention  to  flat  flames, 

then  the  voluminal  stretch  is  constant  and  equal  to  f  .  The  condition  (^k) 

for  b  >  0  can  never  be  realized  but  that  for  b  <  0  can;  if,  when  ,£  is  : 

greater  than  one,  €  is  increased  sufficiently  slowly  for  the  combustion  fiel:' 

2 

to  be  quasi-steady,  then  extinction  occurs  when  it  reachers  -T#/be. 

These  remarks  apply  only  to  the  constant  solutions  in  Fig.  12.  The 
non-constant  solutions  are  available  whatever  the  value  of  €  or  the  sign 
of  b.  No  evidence  was  found  of  quenching  in  the  computations  (cf.  Sec.  2) 
i.e.  the  curves  never  became  tangent  to  streamlines. 
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